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ABSTRACT 

We present the results of the analysis of 7-ray and X-ray data of GRB 050401 taken 
with the Swift satellite, together with a series of ground-based follow-up optical ob- 
servations. The Swift X-ray light curve shows a clear break at about 4900 s after the 
Q ■ GRB. The decay indices before and after the break are consistent with a scenario of 

I continuous injection of radiation from the 'central engine' of the GRB to the fireball. 

C/3 i Alternatively, this behaviour could result if ejecta are released with a range of Lorentz 

' factors with the slower shells catching up the faster at the afterglow shock position. 

The two scenarios are observationally indistinguishable. 
, , The GRB 050401 afterglow is quite bright in the X-ray band, but weak in the 

' optical, with an optical to X-ray flux ratio similar to those of 'dark bursts'. We detect 

a significant amount of absorption in the X-ray spectrum, with Nh = (1-7 ± 0.2) x 
10^^ cm~^ at a redshift of z = 2.9, which is typical of a dense circumburst medium. 
Such high column density implies an unrealistic optical extinction of 30 magnitudes 
if we adopt a Galactic extinction law, which would not be consistent with the optical 
detection of the afterglow. This suggests that the extinction law is different from the 
Galactic one. 

Key words: Gamma-Ray Bursts 



1 INTRODUCTION 

Observations of Gamma-Ray Bursts (GRBs) have shown 
that they are followed by fading X-ray, optical and radio 
afterglows. These are thought to arise when the burst ejecta 
interact with the surrounding medium and produce a shock, 
which propagates in the medium and heats the electrons. 
The latter, cooling by synchrotron emission, produce the 
observed radiation. Studies of afterglows can provide invalu- 
able information on the central engine of GRBs and on the 
circumburst medium, and can possibly distinguish different 
subclasses in the GRBs population. 

Some autho rs (see e.g. f ./azzati et alJ l2002l and ref- 
erences therein; iBerger et alJ 2005r and'' Eamb et alj l2005l 



for recent results) have pointed out the existence of a 
subclass of GRB, whose optical emission is at least ~ 2 
magnitude below that of the average of the optically 
detected bursts. Different models have been proposed to 
explain these "dark bursts", ranging from a cosmological 
origm (iBromm fc L ocb 20021; iFruchteJ Il999^ to scenarios 
where they go off in relatively dense and highly absorbed 
regions (Lazzati et al. 2002). It is also possible that many 
dark bursts may be intrinsically weak so urces, the faint tail 
of the GRBs luminosity distribution jPe Pasq uale et al 
12003 ;} , or sources with a very rapid decav (iGroo^^l 
Il998i) in the optical band. 
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The afterglow emission is seen to decay over time. In 
some bursts, a clear light curve steepening is observed after 
an interval of order days. The break is achromatic and typ- 
ically attributed to the fact that the energy release is coUi- 
mated in a jet. Other irre gular temporal features ar e some- 
times seen in bursts (see IZhang fc MeszarosI l2002l Zhang 
et al. 2005, Burrows 2005 and references therein): exam- 
ples include a rebrightening in GRB 970508, wiggles in 
GRB 020104 and step-like features in GRB 030329. Also, 
"bump" features have been observed in several cases (e.g. 
GRB 970228, GRB 970508, GRB 980326, GRB 000203C) 
and various inter pretations have been p roposed, e.g., "re- 
freshed s hocks" jP anaitcscu ct al. 1998), s upernova com- 
ponents Jbi oom_et_^jl^99i^_jR^idmrtA99^ iGalama et alJ 
I2OOOI). dust echoes llEsin fc BlandfordI l2002^ or microlens- 
mg iGarnavich et alj|2000ir On the other hand, signatures 
in the GRB lightcurve detected at earlier times may provide 
diagnostic information about the nature of the injection and 
eventually probe whether the energy is released impulsively 
during the event or more c ontinuously during the im mediate 
post-burst epoch (see e.g. IZhang fc Meszaro jl2002^ . 

Until recently, most follow-up observations did not start 
until a few hours after the GRB, when the afterglow had al- 
ready faded significantly. This situation has changed with 
the launch of the Swift mission, which provides both a rapid 
alert of GRB triggers to ground-based observers, and rapid 
X-ray and optical/UV follow-up observations of the b urst af- 
terglow. The Swift observatory llGehrels et alJl2'005^ carries 
three science instruments: the Burst Alert Telescope (BAT; 
Barthelmy et al. 2005), which lo cates GRBs with 3' accu- 
racy, the X-ray telescope (XRT, iBurrows et a l. 2005) and 
the U ltra- Violet Optical Telescope fUVQT. [Roming et al^ 
|20053). When BAT detects a GRB trigger. Swift slews to- 
wards the source position within a few tens of seconds. 
Therefore, observations with the Swift instruments yield 
high quality data and cover the poorly investigated epoch 
occurring minutes after the burst. Interestingly, many GRBs 
localized and observed by Swift have shown no optical coun- 
terpart, even when optical observation started 100-200 s af- 
ter the GRB onset. This prov ides evidence for a population 
of "intrinsically" dark GRBs llRoming et al.ll2005htl . 

In this paper, we report the properties of the optically 
faint Swift GRB 050401, and discuss them in the light of the 
current models and scenarios of GRBs. 



2 ANALYSIS OF THE 7-RAY AND X-RAY 
DATA. 

GRB 050401 trig gered the BAT instrument at 14:20:15 UT 
on April 1, 2005 fearbier etaPbOOSl lAngelini et alJl2005ll . 

The refined BAT position is RA=16''31'"16% Dec =2°ll'35" 
with a position uncertainty of 3' (95% C.L., Saka moto et alJ 
[200^. The 7-ray band lightcurve started 9 s before the BAT 
trigger time and it shows 4 main p eaks (see Figure pi. The 
peak count rate was 5000 counts/s iBarbier et al ]|2005). 

Analysis of the BAT data (15-350 keV energy band) 
yields a GRB duration tgo = 33 s. 

We use the mask-weighted technique to subtract the 
background in the BAT for spectral analysis, which is only 
effective up to 150 keV. Swift began to slew towards the 
source about 25 sec after the trigger, while the prompt emis- 



sion was still active. We created separated BAT spectra and 
response matrices for the pre-slew and slew phases, and we 
fitted them jointly (see figure 2) with a simple powerlaw 
model. No significant improvement in is found with a 
cutoff power-law model. Results are reported in the first en- 
try of Table l.\ 

The fiuence detected by BAT in the 15 - 150 keV range 
is (8.6 ± 0.3) X 10~^ erg cm~^. Assuming a redshift of 
z — 2.9 (see later) and a spherically symmetric emission, 
this corresponds to a 7-energy release of 9.6 x lO^'^ erg be- 
tween 15 and 150 keV in the cosmological rest frame of the 
burst (derived by means of the 'k-correction' of Bloom et 
al. 2001). This result differs from that obtained by Chin- 
carini et al. (2005), who reported an energy spectral index 
/3 = -0.13 ± 0.05 (using the 20-150 keV data) and a 7- 
energy release (~ 2.8 x 10^"' erg) in the 15-350 keV band. 
However, we note that GRB 050401 prompt emissio n was 
also detected by Konus- Wind llGolenetskii et al]l2005^ . The 
Konus data (see below) suggest a steepening of the spectral 
index above Eo ~ 150 keV, which might explain the large 
difference between the two analyses. 

As observed by the Konus- Wind instrument, 
GRB 050401 had a duration of ~ 36 s and a fiuence 
of (1.9 3 ± 0.04) X 10"^ erg c m~^ in the 0.02-20 MeV 
band. iGolenetskii et all ll2005l) analyzed the spectrum 
gathered in the first 3 peaks (0-17 s after the trigger, first 
segment) and last peak (24.8-32 s after the trigger, second 
segment) separately (see again Table 1). We note that, if 
we adopt a cutoff powerlaw-model and fix the break energy 
-Eo ~ 150 keV (as inferred from the Konus data), the 
spectral index obtained by Swift, which is averaged over the 
whole gamma-ray emission phase, is consistent with that 
obtained by Konus. 

Observations were not possible with the Swift UVOT 
telescope because of the presence of a 4*'' magnitude star in 
the field of view, while XRT started observations about 130 s 
after the trigger. An unknown bright X-ray source was de- 
tected at R.A. = 16*" 31™ 29^ Dec = 02° 11' 14", within 4 2 
arcseconds of the initial BAT position (^ Angelini et all2005^ : 
the coordinates were later confirmed by ground processing. 
This source subsequently faded, indicating that it was the 
X-ray counterpart of GRB 050401. 

T he first data were taken with the XRT jHill et alJ 
l2005h in its Imaging Mode (IM), then followed by a seg- 
ment of data in the PhotoDiode mode (PD). After that, 
because of the bright star near the edge of the XRT field of 
view, the detector continuously switched between Windowed 
Timing (WT) and Photon Counting (PC) mode. The initial 
PC mode data were piled-up and we did not include them 
in the analysis. In the first day of observation. Swift ob- 
served the X-ray afterglow until 5.6 hours after the trigger. 
More follow-up observations were performed 4.43-6.38 and 
6.48-12.33 days after the trigger. The total exposure time is 
about 48 ks, divided into 21 ks in WT mode and 27 ks of 
PC mode (including piled-up data). 

Analysis of XRT data was performed using the XRT 
pipeline software. The accumulated DN in the IM data were 
converted to a count-rate following the method of Hill ct all 



^ Throughout this paper, we report errors at la, unless otherwise 
indicated 
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I2OO5I . The PD mode data points have been obtained by 
subtracting the contribution of the corner cahbration 
sources. As for the WT mode data, the extraction regions 
for the source and for the background consist of boxes of 
60x40 pixels. The source and background extraction regions 
for the PC data are circles of 20 pixel and 60 pixel radius, 
respectively. We considered the 0-2 grade events for the WT 
mode and 0-12 for the PC mode, which cover events up to 
4 pixels in size. For the spectral analysis, we generated the 
physical ancillary response matrices with the task xrtmkarf, 
while the response matrices (RMs) were retrieved from 
the latest Swift cahbration database, CALDB 20050601 
(http:/ /heasarc. gsfc.nasa.gov/docs/heasarc/caldb/swift/ ). 
At the time of writing, these RMs are considered as 
the most reliable, basing on comparison with spectra of 
calibration sources. We only considered data taken with 
CCD detector temperature T< -48 C. 

We restricted both the spectral and temporal analysis 
to counts within the 0.4-10 keV band since, at the time 
of writing, a good calibration of both WT and PC data 
is available in this band. Data were rebinned in energy by 
requiring a minimum of 15 counts per bin. 

The X-ray lightcurve is shown in Figure |3 The pres- 
ence of a break is clearly evident: if we try to fit the whole 
lightcurve with a single powerlaw, we get an unsatisfactory 
result (xl = 904/114), while the use of a broken powerlaw 
model {At"^ for t ^ tt; At'^^~°'^t°'^ for t ^ U) gives a statis- 
tically significant improvement (xS = 129/112). In this case, 
the best fit parameters are: decay indices qi = —0.63 ±0.02 
and a2 = —1.46 ± 0.07 (before and after the break, respec- 
tively) and break time th = 4900 ± 490 s^. The break time 
is consistent with that reported by Chincarini et al. (2005), 
who report ~ 1300 s corrected for cosmological time dila- 
tion, i.e. their break time has been divided by (1 -\- z). We 
also tried to fit the li ghtcurve with a smooth ly joined bro- 
ken power law model llBeurmann et allll999^ . to determine 
the "sharpness" of the X-ray lightcurve break. However, we 
found that data do not allow us to discriminate between a 
smooth and a sharp transition. 

A fit of the WT spectral data taken before the lightcurve 
break time with an absorbed power-law model reveals a 
considerable amount of absorption, Njf* — (1.5 ± 0.1) x 
lO'^^ cm"'^, clearly in excess of the Galactic value reported 
by Dickey & Lockman (1990), iVg"' = 4.85 x 10^° cm'^ 
and the value inferred from the Galactic extinction map of 
Schlegel et al. (1998), Ng"' = 3.6 x 10^° cm"^ (in this case 
the extinction, Ky = 0.2 has been converted to N^"' using 
results by Predehl and & Schimtt 1995). 

In order to better estimate the column density cor- 
responding to the circumburst medium only, we repeated 
the fit by accounting separately for the Galactic and extra- 
galactic absorption columns, the latter rescaled at z = 2.9. 
We fixed the Galactic absorption to the Dickey & Lockman 
(1990) value. This gives an extragalactic column density of 
Nh = N'lf'^ - N^"' (1.7 ± 0.2) X 10^^ cm-2 (see Table 1; 
spectrum and best fit model are shown in Figure 2J- To as- 
sess the robustness of this detection, we then repeated the 
fit without adding the extragalactic component, obtaining 



^ We note that at the time of the slope change Swift was not 
observing the GRB. 



a xl ~ 405/262. This means that accounting for the extra 
absorption Nh produces a statistically significant improve- 
ment. The F-test has also been widely used to test the sig- 
nificance of a spectral component, although, when applied 
to parameters such as (which is bounde d to be > 0) this 
may be inappropriate in a strict sense fsee iProtassov et alJ 
[2002). For completeness, we report that we checked the F- 
test statistic, obtaining a value of 117 with a probability 
that the improvement is due to chance of ~ lO"^'^. 

We find no evidence for spectral evolution: parame- 
ters consistent with those given above are obtained when 
fitting the spectra taken before and after the 4900 s 
break. An analysis of the PC spectrum with the same spec- 
tral model also does not show any significant difi'erence in 
the spectral parameters (see Table 1). When compared to 
the X-ray afterglows detected by other observatories, like 
BeppoSAX, Chandra and XMM-Newton, the afterglow of 
GRB 050401 appears to be a moderately bright source: the 
1.6-10 keV X-ray fiux normalized at 11 h ours after the burst 
is (2.3 ± 0.2) X 10"^ ^ erg cm'^s"^ fsee lBerger et al.ll200l 
iRoming et ahlbOOSbl De Pasquale et al. 2005). 



3 FOLLOW-UP OBSERVATIONS IN 
DIFFERENT ENERGY BANDS 

Although the Swift UVOT could not observe at the GRB po- 
sition, a series of ground-based optical follow up observations 
was performed, triggered by the prompt Swift localization. 

The afterglow was first identified by Hlvkoff ct <2 
(2005), who detected it in images taken shortly after the 
GRB. Later, following the distribution of the BAT trigger, 
the BAT error circle was imaged with the 40-inch telescope 
at Siding Spring Observatory. These observations consisted 
of two unfiltered 120 s exposures, followed by 41 _R-band 
240 s exposures. Data showed a faint fading source, at co- 
ordinates RA=16:31:28.81, Dec=-|-02:11:14.2 (J2000), i.e. 
within the XRT error circle. This object was not presen t 
on the Digitised Sky Survey plates jMcNaueht et aljlioosll . 

The rapid distribution of the afterglow position, along 
with a finding chart, enabled othe r follow-up observat ions 
including some spectroscopic ones. iFvnbo et alJ (|200^ re- 
ported the detection of a system of absorption lines in the 
optical spectrum detected with F0RS2 at the Very Large 
Telescope, indicating a redshift of z — 2.9. Other obser- 
vations followed, up to ~ 0.5 days after the GRB onset. 
The full list is summarized in Table 2 and the optical light 
curve is shown in Figure |21 We found that the optical decay 
law is not consistent with a single powerlaw (x^ = 24.8 for 
13 d.o.f.). We caution that the optical light curve has been 
obtained by observations performed by difi^erent facilities, 
which used diverse calibrations. In particular, the photom- 
etry of all exposures taken with the Siding Spring Observa- 
tory is relative to the first _R-band image, which was then 
calibrated to match the USNO-A2.0 catalogue red magni- 
tude. Because of this, although the relative magnitudes are 
accurate, the overall normalization can only be considered 
to be accurate at the ~ 0.5 mag level. ARIES measurements 
also refer to USNO-A2.0 stars, while the Maidanak estimate 
is based on stars of the similar USNO-B2.0 catalogue, and 
D'Avanzo et al. 2005, and Greco et al. 2005 data used the 
Landolt catalogue. Therefore some systematic bias may well 
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be present. However, taken at face value the combined data 
suggest a non-monotonic decaying behaviour. We note that 
the decay slope in the first 10000 s (first part of the Siding 
Spring data), is uq — —0.68 ± 0.06, consistent with the first 
X-ray decay. The Siding Spring data are not consistent with 
being constant. A powerlaw fit to these data in isolation 
gives a = 0.56 ± 0.2. The case a = can be rejected at the 
2.5a confidence level. 



4 DISCUSSION. 

4.1 The break in the X-ray hght curve. 

As we have seen, the X-ray light curve of GRB 050401 
registered with XRT shows an initial decay slope of ai — 
—0.63 ± 0.02, which steepens at ~ 4900 s. After this time, 
the source fades with a slope of —1.46 ± 0.07. 

In principle, there are several physical mechanisms that 
can dictate the physical behavior of the afterglow dur- 
ing its initial decline and can produce a break at these 
earl y times. For instance , as described bv lSari et alJ (Il999l) 
and iRhoads et alJ (|l99fll ) , a break in the light curve is ex- 
pected if the fireball outflow is not spherically symmetric 
but coUimated within a jet. Such an effect has been in- 
voked to explain breaks in Gamma-Ray Bur st afterglow light 
curves in several cases (see, for example, | Sari et all Il999l 
IPanaitescu fc Kumaiir200l iFrail et alj|2003h . However, we 
suggest that this explanation is not applicable to the case at 
hand, for at least two reasons. First, according to the fireball 
model of GRB afterglows, the observed temporal decay in- 
dex and spectral slope s hould be linked th rough the so-called 
closure relation (see e.g. IPrice et al.ll2002l) . This relation de- 
pends on the kind of expansion (spherical or jet), on the 
density profile of the medium (uniform or wind) and on the 
cooling state of the electrons responsible for the synchrotron 
emission. In particular, for jet-like emission it should be 
a = 2/3 — 1 or a = 2/3 depending on whether the cooling 
frequency t'e is above or below the X-ray frequency vx (see 
ISari et al .'igg?!. In the case at hand, we find f3 = -0.9±0.03 
(the spectral index is not supposed to change after the begin- 
ning of the jet expansion pheise) and a — — 1.42±0.07 (decay 
index after the break), so none of the above closure relations 
is satisfied. Second, if interpreted within the jet break sce- 
nario the observed parameters are not compatible with the 
relation between the peak energy and the (coUimation cor- 
rected ) 7-ray energy release proposed by iGhirlanda et alJ 
i2004r) . In fact, by using E~,^iso = 3.5 x 10*^ erg (as inferred 
from Konus data) and tt — 4900s we obtain a jet openin g 
angle 6 ~ 1° (see expression 1 in IGhirlanda et al.ll2004) . 
In turn, this corresponds to a coUimation corrected energy 
Ej — 4.8 X 10*^ erg and to a peak energy in the cosmological 
rest frame of Ep^ak = 56 keV, in clear disagreement with the 
value inferred by Konus data (~ 500 keV). 

Another potential reason for a break in the light curve 
is that the energy release is spherically symmetric, but the 
X-ray observation occurs while vc passes through the X- 
ray band, causing the light curve to steepen by 0.25 (see 
fSari et al. 19991. However, in this case we should also ob- 
serve a steepening of 0.5 in the spectral slope while there is 
no evidence for that across the break of GRB 050401. More- 
over, a change of 0.25 is not sufficient to account for the 
steepening in the light curve. 



On the other hand, we find that both the initial shal- 
low decay and the break can be explained by a model in 
which the central engine continues to inject radiative en- 
ergy into the fireball for several t housands of second s . This 
sce nario has been invest i gated bv lZhang fc Meszarod ll200j) 
and lZhang fc Meszarod i2002h . By assuming a source lumi- 
nosity law of the kind L oc , where t is the intrinsic time of 
the central engine (or the observer's time after the cosmo- 
logical time dilation correction), these authors found that 
continuous injection of energy influences the fireball and the 
observed light curves as long as q > — 1. In this case, the 
spectral and decay slopes are linked through the relation: 

a = (l-g/2)/3 + g + l, (1) 

which holds if the observed X-ray frequency is between the 
synchrotron peak frequency and the cooling frequency (see 
later). By using the observed values of a and /3, we obtain 
q ~ —0.5 and g ~ — 1 before and after the break time, re- 
spectively. The latter is consistent with no injection (since 
q < ~1 does not influence the flreball dynamics), while be- 
fore the time break the central engine injects energy with 
a luminosity law L oc t~°'^. The change in the decay slope 
occurs at the point when the central engine ceases to in- 
ject signiflcant amounts of energy. We note that the decay 
slope of the optical and the X-ray flux before the breaks are 
consistent within the errors. This is in agreement with the 
continuous injection model, as long as the optical and the 
X-ray band belong to the same spectral segment. 

There is variant to this scenario, which is observation- 
ally undistinguishable, i.e. a model in which the central en- 
gine activity is as brief as the prompt emission itself but, at 
the end of the prompt phase, the ejecta are released with dif- 
ferent velocities (Lorentz factors, see Panaitescu et al. 2005). 
The fastest shells initiate the forward shock, decelerate, and 
are successively caught by the slowest shells. The consequent 
addition of energy in the blast-wave mitigates the decelera- 
tion and the afterglow decay rate. Assuming that the mass 
M of the ejecta follows the law 

A/(> 7) « 7% (2) 

where 7 is the Lorentz factor, one can flnd an effective s 
value that mimics the effect of non vanishing q index in the 
luminosity law. By following Zhang et al. (2005), this is: 

s = -(10 + 7g)/(2-g), (3) 

therefore, the value q = —0.5 inferred above is equivalent to 
an s-index of s = —2.6. This explanation has been proposed, 
for example, to explain the initial mild decline of the optical 
light curve of GRB 010222, which shows a decay slope of 
ao = -0.7 (Stanek et al. 2001, Bjornsson et al. 2002) for ~ 
10 hours after the trigger, followed by a break and a steeper 
decay. Bjornsson et al. (2004) also proposed that injection of 
energy by slow shells could explain the wiggles in the light 
curve of GRB 021004. 

In the case of GRB 050401, it is noteworthy that a pos- 
sible optical rebrightening seems to take place shortly after 
the change of the slope of the X-ray. In the framework of the 
continuous energy injection model, this could be explained 
by the onset of a "reverse shock". The basic idea is that, 
after the end of the injection phase, a reverse shock wave 
crosses the whole ejecta, heating them and causing a peak 
in the emission. After that, the shocked ejecta start to cool 
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adiabatically once again. However, in order to assess this is- 
sue a more detailed model investigation is required, which 
is beyond the scope of this work. 

As previously discussed, the analysis of the X-ray data 
taken after the break is consistent with a scenario in which 
a "standard" fireball expands in a constant density medium, 
provided that the observed X-ray band lies between the syn- 
chr otron peak frequ ency and the cooling frequency. FoUow- 
ing lSari et alJ lllflflfll ). the closure relation should in this case 
read a = (3/2)/3, which is satisfied within la. One possible 
complication is that, if the fireball expands in a medium 
with constant density, the cooling frequency is expected to 
decrease with time according to uc oc t~^^^. Accordingly, in 
most X-ray afterglows, the cooling frequency is already be- 
tween the peak frequency and the X-ray observing frequency 
less than 1-2 days after the GRB onset. In contrast, our data 
seems to suggest that in the case of GRB 050401 the cool- 
ing frequency remains above the X-ray frequency for about 
10^ s. While this may be explained in terms of relatively low 
values of magnetic field energy (e_g) and density (see Sari et 
al. 1999), we also note that a transit of the cooling frequency 
through the X-ray band after 20,000 s cannot be completely 
excluded by our data: the change in the decay slope and in 
the spectral slope would be 0.25 and 0.5 respectively, and 
hence difficult to detect due to very low statistics in the late 
time XRT detections of the afterglow. 

4.2 The optical and X-ray properties of 
GRB 050401. 

The GRB 050401 afterglow is quite bright in the X-ray band, 
but weak in the optical, with an optical to X-ray fiux ratio 
similar to those of 'dark bursts'. 

In order to compare its properties with those typically ob- 
served in other GRBs, we show in figure |3 the relation be- 
tween the optical and X-ray fluxes for a series of GRB af- 
terglows det ected by BeypoSAX . As poin ted out by severa l 
authors fe.g. lDe Pasauale et alj Eo03. Ro ming et alJl20o'5bl) . 
the GRBs with optical counterparts exhibit a correlation be- 
tween the fluxes in these two spectral bands. Several GRBs 
without optical counterparts also show an X-ray emission 
consistent with that expected by assuming the validity of 
the same correlation, indicating that they may well be sim- 
ply the faint tail of the same population. However, there is 
evidence for a sample of dark GRBs which have "normal" X- 
ray fluxes but with tight upper limits in the optical. Jakobb- 
son et al. (2004) reached similar conclusions by comparing 
the spectral index between the optical and the X-ray band. 
Pox, with the expectations of the fireball model, which re- 
quires Pox ^ —0.5. They found that at least 10 % of the 
events in their sample had Pox > —0.5, and called them the 
truly dark GRBs. To explain the optical faintness of these 
bursts, two main scenarios have been proposed. The first 
idea is that they occur at very high redshift , possibly follow- 
ing th e death of Population II and HI stars iBromm fc Loebl 
I2OO3), in which case the optical flux is washed out by the 
intervening Ly-Q forest. The second idea is that dark GRBs 
have lines of sight passing through large and giant molecu- 
lar clouds (hereafter GMCs). GMCs are rich in dust, which 
extinguishes optical and UV light very efficiently. 

GRB 050401 appears to be an 'intermediate' case. An 
optical afterglow is detected, but is very faint relative to its 



X-ray flux compared to other GRBs with optical counter- 
parts. Its optical to X-ray spectral index is Pox = —0.33, 
which makes this source a dark burst according to Jakob- 
sson et al. (2004) classiflcation. Given the detection of an 
optical counterpart with a likely redshift of 2: = 2.9, we 
can exclude the hypothesis of a very high redshift. On the 
contrary, the X-r ay spectrum indicates a high absorption, 
typical of GMCs llReichart fc Pricdl2002l. iGalama fc Wiierj 
''2001). 

Therefore, a natural question is: can the hypothesis of 
extinction in this medium explain the weakness of the 
optical emission detected? The light detected at Earth 
in the R band (centered at 700 nm) was emitted at a 
wavelength of 180 nm in the GRB cosmological restframe 
at a redshift of 2: = 2.9. The simplest working hypothesis 
is to assume that the same medium is responsible for both 
the X-ray and optical extinction. In this case, the hydro- 
gen column density measured from the X-ray spectrum 
(N_r/ = 1.7 X 10^^ cm~^) would correspond to an absorptio n 
in the V band of Av ~ 10 magnitudes iZombeckl ll99Cf) . 
Assuming a Galactic extinction curve of Ax /Ay oc A~^, 
this results in a predicted extinction of Ax ~ 30 magnitude 
for A = 180 nm. This value is unreasonably high and would 
imply that the optical afterglow was unrealistically bright. 
Thus we can exclude this simple explanation. 

Instead, comparing with the o ptical-to-X-ray flux ratio 
typical of other 'non- dark' GRBs JPe Pasauale et al.ll200l . 
iJakobsson et al .l2004llRoming et alj2005br we could expect 
a plausible extinction of about 3 magnitudes, which is clearly 
not in agreement with the measured X-ray absorption when 
we adopt the Galactic extinction curve. We note that a dis- 
crepancy like thi s has been noted in several other cases (for 
a summary, see IStratta et al.ll20ol) . but for GRB 050401 
we have the advantage of a fairly constrained value of the 
absorption parameter. 

In order to reconcile the value of absorption with the 
likely extinction, a few hypothesis have been proposed. The 
first scenario involves the presence of a gas-to-dust ratio 
much lower than the Galactic one and /or a dust grain siz e 
distribution skewed toward large grains jStratta et alJ2o'o3l . 
This first case is, for instance, typical of dwarf galaxies like 
the SmaU Magellanic Cloud (SMC). In fact, for the SMC 
interstellar medium (ISM) a d ust content ^ 1 /10 of the 
Galactic one has been inferred jPei et alJll992^ . Thus, we 
would have ~ 3 magnitudes of extinction based on the mea- 
sured value of the X-ray absorption. This value is close 
to that required. We note, however, that the best-flt A^j^ 
value has been obtained by assuming Galactic metal abun- 
dance s, while the m etallicity of SMC is 1/8 of the Milky 
Way l|Pei et al.lll99l . If this low metallicity were adopted, 
we would need to adjust the value of Nh upwards by a fac- 
tor of ~ 7, given that the majority of the absorption in the 
X-ray band is produced by heavy elements, so this is not a 
solution. 

Another scenario that could result in a low Av/Nh is 
that there is a change in both the gas-to-dust ratio 
and the distribution of grain size, with the latter en- 
riched in large grains by the effects of the high-energy 
radiation of the GRBs. Dust grains can be heated and 
evaporated by the intense X-ray and UV radiation fields 
up to ^ 20 parsecs from the GRB iWaxman fc Draind 
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l2000l.|Draine fc HaJ2o'oi.lFruchter et al.l20oih . |Perna et'ai] 
ll2003r "and |Pern^^!aaza^ ~ i^2002^ show that the conse- 
quence of exposing dust to intense radiation fields can be 
a grain size distribution flatter than the original one. The 
main reason is that dust destru ction is more efficient on 
small grains. IPerna et al.l (l200.?ll computed the extinction 
curve that is obtained if standard Galactic dust is exposed 
to a GRB and found that the extinction curve can be very 
flat, at least for bursts lasting more than a few tens of sec- 
onds. 

Finally, a distribution of grain size skewed toward large 
grains can also be produced by an efficient mechanism 
of coagulation of smaller grains in a den se environment 
jKim fc Martini Il996l iMaiolino et alJ 120011) . in which case 
the dust-to-gas ratio is unaffected. Both of the above expla- 
nations could lessen the discrepancy between the low UV 
extinction and the high X-ray absorption detected. 



5 CONCLUSIONS. 

We have presented Swift observation of the Gamma-Ray 
Burst 050401 and we have discussed the properties of the 
prompt emission, and, in more detail, the X-ray afterglow. 
The light curve of this burst shows a break 4900 s after 
the trigger, changing from a decay index of ai = —0.63 
to Q2 = —1.46, while the spectral energy index does not 
change. To explain this behavior, we have proposed that the 
'central engine' has been active until the time of the break, 
with luminosity described by the law L oc Another 
possibility is that the central engine activity turned off as 
the prompt emission ended, but the shells had a broad dis- 
tribution of Lorentz factors. In this case, the slowing front of 
the GRB blastwave is continually re-energized by the arrival 
of progressively slower shells and the flux decay is therefore 
mitigated. We find that the decay slope observed before the 
break may be reproduced if the shells were emitted with a 
powerlaw distribution of Lorentz factor, M cc 7~^'^. 
The peculiar behaviour of the optical light curve can be in 
agreement with the two scenarios proposed. 
After the break time, the profiles of the X-ray light curve and 
spectrum are consistent with those expected when a fireball 
expands in a circumburst medium with constant density, 
and the observed X-ray band lies between the synchrotron 
peak frequency and the cooling frequency. 
Even though the X-ray fiux of the GRB 050401 afterglow 
is high, the optical counterpart is faint. This leads to a low 
optical to X-ray fiux ratio similar to that of dark GRBs 
which are likely to be obscured by some mechanism. The 
spectral analysis shows clear evidence of absorption, namely 
Nh = (1.7 ±0.2) X 10^^ cm-^ at the redshift z = 2.9 of the 
GRB. This value is typical of giant molecular clouds where 
star forming regions are located. The detection of a dense 
circumburst medium could lead us to the conclusion that 
this "obscuration mechanism" is, at least in this case, ex- 
tinction. 

Ifowever, the amount of extinction extrapolated by assum- 
ing the Galactic extinction law is far too high to be physi- 
cally acceptable. This may be evidence that the circumburst 
medium is characterized by a dust grain size distribution 
different from the Galactic one, and skewed towards large 
grains. This could be due either to coagulation of smaller 



grains or to small dust grain destruction due to high en- 
ergy photons produced by the GRB. In the latter case the 
dust-to-gas ratio would also be different from the Galactic 
one. 

Acknowledgements: We are grateful to an anonymous 
referee for his/her suggestions that led to a substantial im- 
provement of the draft. APB, KPA acknowledge support for 
this work at Leicester by PPARC. SZ also thanks PPARC 
for support through an Advanced Fellowships. 



REFERENCES 

Angelini L. et al., 2005, GCN 3161 
Barbier et al., 2005, GCN 3162 

Barthelmy S.D. et al., 2005, Sp. Science Rev, in press 
Beuermann K. et al. 1999, A&A, 352, L26 

Berger E. et al., 2003, ApJ, 590, 379. 

Berger E. et al., 2005, ApJsubmitted, l |astro-ph /0505107| 
Bjornsson G., Hjorth J., Pedersen K. et al., 2002, ApJ, 
579L, 62. 

Bjornsson G., Gudmundsson E.H. & Johannesson G., 2004, 

ApJL, 615, 77. 
Bloom J.S. et al., 1999, Nature, 401, 453 
Bloom J.S, FraU D.A., Kulkarni S.R. et al., 2001, ApJ, 508, 

L21 

Bromm V. & Loeb A., 2002, ApJ, 575, 111 
Burrows D. et al., 2005, Space Science Review, in press 
Chincarini et al. 2005, submitted to ApJ, astro-ph/050645^ 
D'Avanzo P. et al., GCN 3171 

De Pasquale M., Piro L., Perna R. et al. 2003, ApJ, 592, 
1018 

Dickey J.M. & Lockman F.J., 1990, ARAA, 28, 215 

Draine B.T. & Hao L., 2002, ApJ, 569, 780 

Esin, A. A. & Blandford, R., 2000, ApJL, 534, L51 

Frail D. 2003, lAU Colloquium 192, Supernovae: 10 Years 

of 1993 J Valencia, astro-ph/0309557| 
Frontera F. et al., 2000, ApJ, 127 59 
Fruchter A.S. et al., 2001, ApJ, 563 597 
Fynbo J. et al., 2005, GCN 3176 
Fruchter A.S., 1999, ApJ, 512L 1 
Galama, T.J., et al., 2000, ApJ, 536, 185 
Galama, T.J. & Wijers, R.A.M.J. 2001, ApJL, 549, L209 
Garnavich, P., Loeb, A. & Stanek, K., 2000, ApJL, 544, 11 
Gehrels N. et al., 2005, ApJ, 621, 558 
Ghirlanda G., Ghisellini G. & Lazzati D., 2004, 616, 331. 
Golenetsku S. et al., 2005, GCN 3179 
Greco G. et al., 2005, GCN 3319 

Groot P.J, Galama T.J., Vreeswijk P.M. et al. 1998, ApJ, 
502L 123 

Hill J.E. et al., 2005, in preparation 

Jakobsson P. et al., 2004, ApJL, 617, L21 

Kim S.H. & Martin P.G., 1996, ApJ, 462, 296 

Lamb, D.Q. et al., 2005, New Astronomy Reviews, 48, 423 

Lazzati D. et al. 2002, MNRAS 

Maiolino R. et al., 2001, A&A, 365, 37 

Misra K. et al. 2005, GCN 3175 

McNaught R. et al., 2005, GCN 3163 

Panaitescu A., Meszaros, P. & Rees, M.J., 1998, ApJ, 503, 
314 

Panaitescu A. & Kumar P., 2002, ApJ, 779, 789 



Swift and ground-based observations of GRB 050401. 



Panaitescu A. et al., 2005, MNRAS submitted, 

astro-ph/0508340 
Pei Y.C. et al., 1992, ApJ, 395, 130 
Perna R. & Lazzati D., 2002, ApJ, 580, 261 
Perna R. et al., 2003, ApJ, 585, 775 
Predehl P. & Schmitt J.H.M.M., 1995, A&A, 293, 889 
Price P. et al., 2002, ApJL, 572, L51 
Price P. et al., 2005, GCN 3164 

Protassov R., van Dyk D.A., Connors A. et al. 2002, ApJ, 
571, 545 

Reichart E., 1999, ApJL, 521, Llll 
Reichart E. & Price P., 2002, ApJ, 565, 174 
Rhoads J. et al., 1999, ApJ, 525, 737 
Roming P. et al., 2005a, Space Science Review, in press 
Roming P. et al., 2005b, to be submitted to Science, as- 
tro/ph0509273 

Rykoff E. S. et al., 2005, accepted by ApJL, 

astro-ph/0508495 
Sari R. et al., 1999, ApJL, 519, 17 
Sakamoto T., 2005, private communication. 
Schlegel D.J., Finkbeiner D.P & Davis M, 1998, ApJ, 500, 

525S 

Soderberg A.M. et al., 2005, GCN 3187 
Stratta G. et al., 2004, ApJ, 608, 846 
Stanek K., Garnavich P.M, Saurabh J., et al., 2001, ApJ, 
563, 592. 

Waxman E. & Draine B.T., 2000, ApJ, 537, 796 
Zahharov. et al., 2005, GCN 3174 

Zhang B., Fan Y.Z., Dyks J. et al., 2005, submitted to ApJ, 

astro-ph/0508321 
Zhang B. & Meszaros P., 2001, ApJ, 552L, 35 
Zhang B. & Meszaros P., 2002, ApJ, 566, 712 
Zombeck M., 1990, Handbook of Astrophysics, Ed. Cam- 
bridge University Press 



8 



Instrument 



Section 



Eq 



(1022) 



Swift: 



BAT Prompt Emission 
WT pre-break 
WT post-breaJj 
PC 



-0.53 ±0.07 
-0.9 ±0.03 
-0.9 ±0.13 
-0.75 ±0.15 



136/115 
1.7 ±0.2 280/261 
1.7 16.3/21 
1.7 21.8/18 



Konus-Wind; 



first segment 
second segment 



-0.15 ±0.16 
+0.17 ±0.21 



-1.65 ±0.31 156 ±45 
-1.37 ±0.14 119 ±26 



Table 1. Values of the parameters for the spectral fit of GRB 050401. Swift data have been fitted with a simple power-law model (energy 
index (3) and divided into four sections; prompt emission, WT data before the break in the light curve, WT data after the break, PC 
data. For the latest two sections, N[{ has been kept fixed to the value obtained from the WT pre-break data. Errors are at 68% confidence 
level. For comparison, we also reported the results of the fitting of Konus-Wind prompt emission data performed by Golenetskii et al. 
(2005). A Band model has been adopted by these authors, so in this case /3, /3i axe the low and high energy indices and Eq is the breaJj 
energy. 



Time after GRB (days) 


R magnitude 


Reference 


0.040 


21.05 ±0.3 


Siding Spring Observatory 


0.043 


21.25 ±0.3 




0.051 


20.85 ± 0.2 




0.055 


21.35 ± 0.3 




0.063 


21.85 ± 0.4 




0.071 


21.55 ± 0.3 




0.086 


21.55 ±0.2 




0.102 


21.85 ±0.3 




0.126 


21.85 ±0.3 




0.149 


21.85 ±0.4 




0.173 


22.05 ±0.3 




0.24 


21.27 ±0.2 


Misra et al. 2005 (GCN 3175) 


0.39 


21.97 ±0.2 


Kahharov et al. 2005 (GCN 3174) 


0.46 


22.47 ±0.4 


Greco et al. 2005 (GCN 3319) 


0.47 


22.95 ±0.1 


D'Avanzo et al. 2005 (GCN 3171) 



Table 2. Log of GRB 050401 Optical observations. Values quoted in this table have been corrected for Galactic extinction. As explained 
in the text, the zero-points for these observation might be uncertain for ~ 0.5 magnitude due to USNO-A2.0 calibration. 
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Maskweighted 1 s Lightcurve 
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15-25 keV 



20 
MET Time - 134058014.784 

Figure 1. GRB 050401 light curve in Gamma-Rays. On the Y-axis we show the background subtracted count/s per fully illuminated 
detector for an equivalent on-axis source. 
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data and folded model 




channel energy (keV) 



Figure 2. GRB 050401 7-ray spectrum detected by BAT before (red colour) and during the slew (black). Both spectra are consistent 
with the same fitting power law model, plotted as a solid line (see first entry in Table 1). 
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Figure 3. GRB 050401 afterglow light curves. The crosses represent the X-ray light curve registered with XRT (0.4-10 keV band). The 
first point is taken in Image mode, followed by 2 Photodiodc points. After that, data have been taken in Windowed Timing (WT) mode 
from 0.14 up to 9 ks after the trigger (section between the two vertical dashed lines), and in Photon Counting (PC) mode from 13.6 up 
to 1050 ks after the trigger. The optical afterglow light curve is plotted in the upper part of the figure. Triangles: Siding Spring data. 
Circle; ARIES (Misra et al. 2005) data. Filled circle: MAO (Kahharov et al. 2005). Square: Bologna (Greco et al. 2005) data. Filled 
square: TNG (D'Avanzo et al. 2005) data. 
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GRB 050401 WT spectrum 
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Figure 4. GRB 050401 X-ray spectrum registered with XRT in WT mode before the breaJf. The solid line is the best fitting absorbed 
power-law model (see text and Table 1 for details). 
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Figure 5. Optical vs. X-ray flux of GRB 050401 afterglow compajred to BeppoSAX GRBs. Open Circles: GRBs with optical counterpart. 
Solid arrows: GRBs without optical counterpart. Dashed arrows: doubtful afterglows. Dotted arrows: upper limits. Short-long dashed 
line: best fit of optical vs. X-ray flux (adapted from De Pasquale et al. 2003). 



